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Abstract 

We present a new approach to the problem of G fe,i -constrained (fc, l < 3) multi-degree reduction of Bezier 
curves with respect to the least squares norm. First, to minimize the least squares error, we consider two 
methods of determining the values of geometric continuity parameters. One of them is based on quadratic 
and nonlinear programming, while the other uses some simplifying assumptions and solves a system of linear 
equations. Next, for prescribed values of these parameters, we obtain control points of the multi-degree 
reduced curve, using the properties of constrained dual Bernstein basis polynomials. Assuming that the 
input and output curves are of degree n and m, respectively, we determine these points with the complexity 
O(mn), which is significantly less than the cost of other known methods. Finally, we give several examples 
to demonstrate the effectiveness of our algorithms. 

Keywords: Constrained dual Bernstein basis, Bezier curves, Multi-degree reduction, Geometric continuity, 
Quadratic programming, Nonlinear programming. 


1. Introduction 

Let 11“ denote the space of all parametric polynomials in of degree at most n; II* := II„. 

A Bezier curve P n £ 11^ of degree n £ N is the following parametric curve: 

n 

P n (t) -Y^PiBiif) (0 < t < 1), (1.1) 

i—0 

where po,Pi, ■ ■ ■ ,p n £ R d are so-called control points , and Bq , B™,..., B ™ are the Bernstein polynomials of 
degree n given by 

:= Qt*(l-*)"-< (0 < i < n). (1.2) 

In this paper, we consider the following problem. 

Problem 1.1. [G fc,, -constrained multi-degree reduction] 

For a given Bezier curve P n of degree n, find a Bezier curve R m of lower degree m, 

m 

Rm(t ) := X G#r(t) (0 <t< 1), (1.3) 

i=0 

so that the following conditions are satisfied: 

(i) P n and R m are G k ’ 1 -continuous (—1 < k,l < 3 and k +1 < m— 1) at the endpoints, i.e., 

^-R m {t) = ^-P n (ip(t)) (t = 0; i = 0,1,..., k), 

cF d j ( 1 ‘ 4 ) 

— R m (t) = — P n (<p(t)) (t = 1; j = 0,1,.... Z), 

where (p : [0,1] —»• [0,1] is a strictly increasing function with (p( 0) = 0 and y?(l) = 1; 
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(ii) value of the squared L 2 ~error 


I p _ j? 

\ J n 1L n 


|2 
11/2 


:= [ (1 - t) a t p \\P n (t) - R m (t)\\ 2 dt (a, /3 > - 1) 

J o 


(1.5) 


is minimized in the space 11^, where || • || is the Euclidean vector norm. 

Problems of the above type have been recently discussed in several papers BHEEEEE, usually un¬ 
der simplifying assumptions <//(0) = ^/(l) = 1, which implied, for example, the hybrid C' 1,1 /G 2 ’ 2 -constrained 
degree reduction, meaning that we impose constraints of G 1,1 -continuity, followed by G 2,2 -continuity, at the 
endpoints. Most of the known algorithms solve a system of normal equations to get control points of the 
multi-degree reduced curve m- Consequently, solution depends on the inverse of a certain matrix, so the 
obtained formulas are not truly explicit and the cost of the method is high (see, e.g., HE E ). For extensive 
lists of references, see the recent papers of Lu [fjj, or Rababah and Mann [E|- The conventional problem of 
degree reduction differs from Problem 11.11 in considering, instead of condition (i), the G fc,z -continuity at the 
endpoints of curves, i.e., 

R$(0) =P}?\0) (t = 0 , 1 ,..., fc), 1 
= (j = 0,l,...,/)■ J 

In the past 30 years, many papers dealing with this problem have been published (see, e.g., BSE EE). 
In particular, in [151] . two of us have proposed a method based on the use of the so-called dual Bernstein 
polynomials , which has complexity O(mn), the least among the existing algorithms. In the present paper, 
we apply an extended version of this method as an essential part of the algorithms of solving Problem 11.11 
Such an approach allows us to avoid matrix inversion. Assuming that — 1 < fc,l < 3 and including the hybrid 
cases, there are 37 continuity cases which require computation of the continuity parameters. Those variants 
of the problem differ, and we have not proven that in each case a unique solution exists. 

The outline of the paper is as follows. Section [2] contains a preliminary material. In Section [3j we relate 
the G fc ’ Z -continuity conditions with the control points of the curves P n and R m . Section [3] brings complete 
solutions of Problem o with and without the simplifying assumptions. Section [5] deals with algorithmic 
implementation of the proposed methods. In Section [Gl we give some examples showing efficiency of our 
methods. Conclusions are given in Section [7] 


2. Preliminaries 


In this section, we introduce necessary definitions and notation. 

We define the inner product (-,-)a:,/3 by 

( f,9)a,p~ [ (1 -t) a t p f{t)g(t)dt (a, j8>-1). (2.1) 

Jo 

There is a unique dual Bernstein polynomial basis of degree n 


D n D n 
u o > u \ > • 


n n g. TT 


associated with the basis & so that 




( i,j = 0, l,...,n), 


where Sij equals 1 if i = j, and 0 otherwise. 

Given the integers fc, l such that fc, l > — 1 and k + l < n — 1, let 11^’^ be the space of all polynomials of 
degree at most n, whose derivatives of orders 0,1,..., fc at t = 0, as well as derivatives of orders 0,1,..., l at 
t = 1, vanish. We use the convention that derivative of order 0 of a function is the function itself. Clearly, 
dim II n’ 1 ^ = n — k — l — 1, and the Bernstein polynomials {B^ +1 , B ^ +2 ,..., } form a basis of this 

space. There is a unique dual constrained Bernstein polynomial basis of degree n 


’-^/c+2 ’* 


y-v(n,/c,Z) 



satisfying the relation / nG.M), bv\ = ( i,j = k + 1, fc + 2,..., n — l — 1). Obviously, we have 

\ J / a,/3 

11 = Gf, which corresponds to the case without any constraints. For properties of the polynomials 
Bf and G", see [E, Appendix A]. 
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Forward difference operator is given by 

A ° qi :=qi, A fc ® := - A*- 1 ^ (k = 1,2,...). 

We use C p,q /G k ’ 1 notation to describe the hybrid constraints, where p,q £ {—, 1} and (fc > 2 or l > 2). 
In the case of k > 2 and p = 1, we set <p'(0) := 1. Similarly, for l > 2 and g = 1, we set (p'(l) := 1. Setting 
p := —, g := — means that we do not fix <p'(0), <//(!), respectively. Clearly, C~’~/G k ’ 1 denotes G fe,i . 


3. Geometric continuity 

In this section, we relate the G fc,i -continuity conditions OD with the control points. We limit ourselves 
to k, l < 3 cases, which are the most important from a practical point of view. 

Remark that the control points r ±, r ^,..., rk depend on the parameters 

\j ■— ipW ( 0 ) (j = 1 , 2 ,..., k), 

while the points r m _i,r m _ 2 ,..., r m _/ depend on 

Pj := { 1) (j = 1,2, ...,0- 

Now, let us recall the well known formulas (see [?J, also SEU). When k = 3, we have: 

Tl 

ro=Po, rx=p 0 -\ -AiApo, (3.1) 

m 


n 

r2=Po~\ - 


m 


2Ai + 


1 

m — 



Ap 0 + 


(^- 1)2 
(m - 1)2 


A?A 2 p 0 , 


(3.2) 


p 3 = PO H- 

m 


3Ai H- 7 A 2 + 


1 


m — 1 (m — 2)2 


An 


+ 3 


( n ~ 1)2 


Af + 


1 


(m — 1)2 to — 2 


-A 1 A 2 


Ap 0 

fa - 2)3 


a 2 , TO ~ TO 3 ,3.3 

A po + — A : A po. 

(to - 2)3 


(3.3) 


In the case of k = 2, we use (ELD and ( 22 ). For k = 1, formulas (E3D hold. Analogously, when l = 3, we 
have: 


T m —Vni T m— 1 —Pn P'1 ^Pn—lt 

m 


T m —2 — Pn 

m 


2pi — 


TO — 1 


M 2 


. fa — 1)2 o . 2 

Apn-l + 7-7—MlA 2 p„_2, 

(to - 1)2 


(3.4) 

(3.5) 


Tm —3 — Pn 

m 


3pi - 


1 


TO — 1 


+ 3 


fa - 1)2 

(to - 1)2 


Mi 


M 2 

1 


m — 2 


(to - 2)2 
Mi M 2 


M3 


Ap n _i 

fa - 2) 3 


A V-2 - r ;; 3 MiA 3 p„- 3 . 

(to - 2)3 


(3.6) 


In the case of l = 2, we use (E3D and (13.51) . For l = 1, formulas (13.41) hold. 


4. G fc,I -constrained multi-degree reduction problem 

4-1. Multi-degree reduction of Bezier curves with prescribed boundary control points 
First, we discuss the following model problem of constrained multi-degree reduction: 

Problem 4.1. [Multi-degree reduction with prescribed boundary control points] 

Given a Bezier curve P n £ II((, 

n 

P n (t) 

2=0 
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we look for a Bezier curve R m £ 11^ (m < n), 


m 

R m (t):=J2riBr(t), (4.1) 

2=0 

having the prescribed control points ro, 7 * 1 ,..., r/~ and r m _/, , r m; that gives minimum value of the 

error 

:=\\P n -R m \\l 2 = [ (l-t) a t^\\P n (t)-R m (t)\\ 2 dt (a, 0 > —1). (4.2) 

Jo 

Given the points p t := (pn,pa ■ ■ ■,Pid ) £ (i = 0,1,..., n) and r t := (rii,r i2 ... ,r id ) £ R d (i = 0,1,... ,m), 
we use notation p ft , r h for the vectors of hth coordinates of the points po,Pi, ■ ■ ■ ,Pn and ro,ri,... ,r TO , 
respectively: 

P •— \Poh lPihi • • ■ • Pnfa] , r .— [ro/i , rih ,..., r m ^] (6 1, 2,..., g?) . 

As an extension of the result given in [lj| (see also ®), we obtain the following theorem. 


Theorem 4.2. The inner control points r* = ( rn,rt 2 ,..., r^) (fc + 1 < i < m — l — 1) of the curve (14.11) . 
being the solution of the Problem are given by 


ri = (i = k + 1, k + 2,..., m — l — 1), 

3=0 


where 


and 


_/ Tjn i 

VlJ ■— \ n j 7 


a,/3 


-v=«- t ")"fy> u o=o,i,...,«). 


\h =0 h=m—l. 

The squared L 2 - error (S3 is given by 


j-hj\h 


(4.3) 


(4.4) 


(4.5) 


E (a 'P) = ^ {!nn(p h ,p h ) + Imm(r h ,r h ) - 2I nm (p h , r h )} , 


(4.6) 


h -1 


where for a := [ao, ai,..., a at] and b := [bo, &i,..., 6 m], we define 


Inm(&7 b) := 


where B(ct,/3) := is the beta function. 

Proof. Let us write 


B(a + 1, /3 + 1) 
(a + /3 + 2 )at+m 


ee(T)(t) (a + l)jV+M-i-j(^ + l)i+jRi^i, 

i =0 j =0 v v J ' 


where 


Rmit) = S m {t) + T m (t), 


S m (t):= J2 ^(t), T m (t):=[J2+ E 


i=k -\-1 


v, 2=0 i—m—l/ 


Using the degree elevation formula (see, e.g., 0, §6.10]; we adopt the usual convention that (“) = 0 if v < 0 
or v > u) 


7 £ 


h=0 


n — m\ I n 

h-i )\h 


-l 


W), 


we write 


T m (t) = Y / d J B?(t), 

3=0 
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where 


Now, we observe that 




I p _ p 

I J n, 1 


= \\w n - s m \\i = jr l\i - t ) a tP [wZ(t) - stitfdt, 

h=i Jo 


where 


W n (t) := 

Sm(t ) := , 


2—0 


with 


Vi := Pi - di. 


Thus, we are looking for the best least squares approximation for W% (h = 1,2 ,d) in the space 11™’^. 
Remembering that B™ and D^ rn,k ’ l ' > (k + 1 < i < m — l — 1) are the dual bases in the space II m ,l \ we obtain 


n n 

Tj = Y, V :i (b™, = y ^Vj(f>ij (i = k + l,k + 2, ... , m - l - 1), 

3=0 01,13 j =0 


which is the formula (14.31) . 

Proof of (14.61) uses an argument similar to the one given in 


m- 


□ 


Remark 4.3. Let us define the quantities ipij (i = k + 1, k + 2,..., m — l — 1; j = 0, 1, ..., n), related to the 
coefficients <pij (cf. (14.41) 1 by the following formula: 


fm — k — l- 2\ /m\ 1 /n\ (a + l + 2 ) n _j (/? + k + 2 )j 
V i — k — 1 )\i) \j) {& + l + 2)i + i(/3 + k + 2)k+i 3 


(4.7) 


Observe that the quantities ipij can be put in a rectangular table and the entries of this ^-table can be 
computed using Algorithm 4.2], assuming that C\ := k + 1, C2 := l + 1, oq := a and a 2 := /?■ Note that 
the complexity of this algorithm is 0(mn). 


4-2. G k l -constrained multi-degree reduction 

Coming back to the problem of G fc,z -constrained multi-degree reduction (see Problem O). let us notice 
that the formulas m dMt with fixed parameters A; and pj (cf. (J3J constitute constraints of the form 
demanded in Problem 14.II As a result, the control points (14.31) depend on these parameters. 

Now, the optimum values of the parameters can be obtained by minimizing the error function (14.61) . 

E&’ft = £?(“■« (A 1} A 2 ,..., A fe , (4.8) 

depending on {A^} and {pj} via formulas (13.11) (13.61) and (14.31) . 

For a minimum of function (14.81) . it is necessary that the derivatives of with respect to the param¬ 

eters are zero, which yields the system 


E E [F mj (* h ) - Fnj(p h )\ 

h =1 j—u 
d m—v 

E E [Fmj{v h ) - F nj { p*)] 

h—1j=k +1 


° r :> h = 0 
d\ u 

(u = 1,2,. 

..,k), 

^ 3- 

II 

O 

(v = 1,2,.. 

■J), 


where we use notation 


-ftj(q) 


1 

(q + f3 + vn + 



(p H - l)t+m— i—j (/^ “1“ Qi 


(4.9) 
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with q = [q 0l qi,... ,q t }- 

In the case of k = l = 3, we compute the partial derivatives of hth coordinates of the control points m 
ED. We obtain: 


dr. 


ih 


d\i 


= 


% A Poh 

2SAp m + 2A 1 ^l|a 2 p (Jk 


3% A Poh + 
0 


2A 1 + 7 ^ 2 ^] 3 (m-lh A2 P0h + 3A? ("JJ) 3 , A 3 p, 


( m - 2)3 


Voh 


dr ih 

dX 2 


(1 


T=Ij4 A Poh 




Oh 


0 


(* = 2 ), 

(i = 3), 

(i = 0,1; m — 3 < i < m), 


dr. 


ih 


8 X 3 


dr. 


ih 


dpi 


= 


(m-2) 3 A P° h (* — 3), 

0 (* = 0,1,2; m — 3 < * < m), 

~^ A Pn-l,h 

~^m A P n - l > h + 2^1 ( m _!) 2 2 A“ Pn-2,h 

3 m A P n — l,h T 2/Xi — ' A Z 2 P2 3 A 2 Pn-2,h ~ 3p\ ( m -2) 3 A3 Pn-3,h 


dr ih 
dp 2 


> n 

’ (m —1)2 


3 A Pn—l,h ~ 3^1 ' fa _ 2 Vi A 2 p ra - 2 ,ft (i = m - 3) 


dr. 


ih 


dp3 


( 771 — 2)3 

0 

(m— 2) 3 A Pn—l,h (l = m — 3), 

0 (0<i<3;m — 2<i< m). 


(i = m — 2), 

{i = m — 3), 

(0 < i < 3; i = m — 1, m), 


(i = 1). 

(i = 2 ), 

(* = 3), 

(* = 0; m — 3 < i < m), 

(4.10) 

(4.11) 

(4.12) 

(i = to — 1 ), 

(i = m — 2 ), 

(i = m — 3), 

(0 < i < 3; * = m), 

(4.13) 

(4.14) 

(4.15) 


Notice that the partial derivatives of /ith coordinates of control points (14.31) depend on (14.101) (14.151) in 
the following way: 


dr ih 

dX u 


E 

3=0 


n — m 


-1 k 




(4.16) 


dr ih 

dp v 




-1 m—v 

E 

g—m—l 


n — m 
j ~ 9 


b-.foj !i 

a#*,, ■ 


(4.17) 


One can easily see, that when k, l < 3, we compute by (14.161) . (14.171) if k < i < m — l, and 

by (14.101) (14.151) otherwise. Finally, we put the expressions (14.101) (14.171) into the equations of system (14.91) . 

Observe that for k > 2 or l > 2, system (HU) is nonlinear, which makes it quite difficult to solve. 
Furthermore, from a practical point of view, we additionally require that X\,pi > 0, which results in the 
same directions of tangent vectors at the endpoints of curves ED and ED- Therefore, to guarantee that 
these conditions will be satisfied, it is not enough just to solve the system ED- 

Now, let us discuss two possible ways of determining the values of geometric continuity parameters. 


4-2.1. Determining the G k ' 1 parameters using optimization methods 

It is easy to check that if (k = 1 and l < 1) or (2 = 1 and k < 1), then the error (14.61) is a quadratic 
function of continuity parameters. 
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In the case of (k = 2 and l < 2) or (7 = 2 and k < 2), the error (14.611 is a fourth-degree polynomial 
function of continuity parameters. 

For (k = 3 and l < 3) or (7 = 3 and k < 3), the error (14.61) is a sixth-degree polynomial function of 
continuity parameters. 

To find the optimum values of parameters Xi,jm in the case of G 1,1 -constrained multi-degree reduction 
problem, assuming that a, j3 = 0, Lu and Wang [l0( solve the quadratic programming problem, subject to 
the constraints 

Ai > do, p\ > d\, (4-18) 

where do and d\ are positive lower bounds, usually prescribed to small values (they set 10 -4 for both lower 
bounds in the examples section). Such approach can be used in the cases which result in the quadratic error 
function m - One can solve the quadratic programming problem using, e.g., an iterative active-set method, 
which is implemented in many software libraries. The active-set mechanism used by standard quadratic 
solvers is described in 0 §6.5], 

Analogously, one can observe that for k = 2,3 or l = 2,3, the problem of minimizing the error (14.61) . 
subject to the constraints (14.181) is a nonlinear programming problem. To solve it, one can use, for instance, 
a sequential quadratic programming (SQP) method (see, e.g., [l|, §15-1]), which is available in many software 
libraries. 

4-2.2. Determining the C p,q /G k ' 1 parameters by solving a system of linear equations 

In the case of G 2,2 , Rababah and Mann [TjJ simplified the problem by considering G 1,1 -continuity at the 
endpoints, i.e., they set Ai = pi := 1. Later, this approach was also used by Lu @. In fl3], the same idea 
was used to simplify the G 3,3 case, and the authors noted that such approach leads to a system of linear 
equations. 

Now, we generalize the above-described approach for any k. I such that — 1 < k, l < 3. If k > 2, we 
set Ai := 1, which implies G 1 -continuity at t = 0 and consequently, G k ’ 1 constraints become C 1 ’ q /G k ’ 1 
constraints, where q £ {—, 1}. Similarly, when l > 2, we set pi := 1, which implies G^-continuity at t = 1 
and consequently, G k ’ 1 constraints become G p, 1 /G fc,z constraints, where p £ {—, 1 }. 

Notice that in the cases of k = 2, 3 or l = 2, 3, the above-described method leads to the linear system (Oil 
and the error (imp is a quadratic function of the continuity parameters. However, in the cases of k = 1 or 
1 = 1, there is no guarantee that the solution satisfies Ai > 0 or pi > 0, respectively. In the case of solution 
with nonpositive values of these parameters, our choice is to solve a quadratic programming problem, subject 
to the constraints with prescribed positive lower bounds for the parameters (see (14.181) 1. 

Observe that this approach uses no simplifying assumptions for k, l < 1. 

Let us denote the above-described approach to Problem ll.ll as C p,q /G k ' 1 -constrained multi-degree reduction 
of Bezier curves. 

5. Algorithms 

In this section, we show the details of implementation of the proposed method of G fc,, -constrained multi¬ 
degree reduction of Bezier curves. Moreover, we give a short description of G^/G^-constrained multi-degree 
reduction algorithm. 

5.1. G k ’ 1 -constrained multi-degree reduction algorithm 

Now, we give the method of solving Problem ll.il summarized in the following two-phase algorithm. 

Phase A of the algorithm consists in finding values of the parameters Ai and pj to minimize the error (14.61) . 
which—by the results given in Theorem l4.2l — depends only on these parameters. The idea is based on solving 
the quadratic or nonlinear programming problem (see 94.2.11) . Notice that when k, l < 1, we can compute the 
(fij coefficients (cf. (14.71) 1 and omit the remaining steps of Phase A, since there are no continuity parameters 
to determine. During Phase B, we use the results of Theorem 14.21 and the obtained values of continuity 
parameters to compute control points ro,r\,..., r m . Most of the known algorithms solve a system of normal 
equations, to get the inner control points of multi-degree reduced curve m- Such approach makes these 
points dependent on the inverse of a certain matrix. Our formulas do not require matrix inversion. What is 
more, the complexity of Phase B is 0(mn), which is significantly less than the cost of other known methods 
for this phase. The algorithm works for any k and l not exceeding 3. 

Algorithm 5.1. ]G k ’ 1 -constrained multi-degree reduction] 

Data: a, /3 - parameters of the inner product ED; 
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n, PchPi, ■ ■ ■ iVn - degree and the control points of the Bezier curve (HD; 
m - degree of the reduced Bezier curve OD; 

k, l - orders of the G-continuity at the endpoints of the curve OD; 
do, d± - lower bounds for the parameters Ai and p\, respectively (cf. 

Assumptions: n > m > 0; —1 < k,l < 3; k + l < m — 1; do, d± > 0; a, /3 > —1. 

Result: control points ro,ri ,..., r m of the G kl -constrained multi-degree reduced Bezier curve OD- 

Phase A 

Step I Compute {4>ij} (i = k+1, k+2,..., m — l — 1; j = 0,1,..., n) by fs. Algorithm f.2] and formula (14.71) 
(see Remark \4.3\ ). 

Step II Check if the remaining steps of Phase A can be omitted: 

If (k < 1) and (l < 1) then go to Step VI. 

Step III Compute E^ a 'd\Xi, \ 2 ,..., ■ ■ ■, Pi) by (14.61) . 

Step IV Determine set c of constraints: 
c := {Ai > d 0 , /ii > di}; 

If (k < 1) then c:= c \ {Ai > do}; 

If (l < 1) then c := c \ {pi > di}. 

Step V 

If (k > 1 or l > 1) then 

obtain Ai, A 2 ,..., A k, emd /ii, P 2 , ■ ■ ■, m by solving the nonlinear programming problem of minimiz¬ 
ing the error OD> subject to the constraints c; 

else 

obtain X\, X 2 , ■■■, Xk, and p±, p. 2 , ■ ■ ■, Pi by solving the quadratic programming problem of minimiz¬ 
ing the error OD, subject to the constraints c. 

Phase B 
Step VI Compute 

1. r 0 ,n,...,r k by (HTTP (OD ; 

2. r m -i,r m -i + i,...,r m by OD (OD . 

Step VII Compute Vo,Vi,... ,v n by OD. 

Step VIII Compute r k+ i,r k + 2 , ■■■ ,r m -i-i by (OD . 

Step IX Return the solution, i.e., the control points J"o,n, ■ ■ • of the reduced Bezier curve OD- 

5.2. C p ’ q /G k ’ 1 -constrained multi-degree reduction algorithm 

Now, let us give the outline of the two-phase C' p,9 /G fc ’ i -constrained multi-degree reduction algorithm. 
Phase A of the algorithm implements the ideas discussed in M.2.21 therefore, it solves the system of linear 
equations (14.911 to determine values of the continuity parameters. In the case of solution with nonpositive 
values of Ai or p\, which can happen when k = 1 or l = 1, the algorithm solves a quadratic programming 
problem, subject to the constraints with prescribed positive lower bounds for the parameters (see ( 4.181) 1. 
An example of a resulting Bezier curve that does not satisfy the positive condition can be found in 6j, Fig. 
1(a)]. We performed more than 40 different tests (results of some of them are available in the next section). 
None of them caused such problem. Phase B is the same as for Algorithm 15.11 The algorithm works for 
any k and l not exceeding 3. For details, see our implementation in Maple™ 13 available on the website 
http://www.ii.uni.wroc.pl/~pgo/GDegRed.mws. 

Obviously, Algorithm 15.11 costs more, but also produces more accurate results, since for the C p ’ q /G k,l ‘- 
constrained approach we additionally assume that Ai = 1 when k > 1, and fi\ = 1 when l > 1. 








6. Examples 

This section provides of the application of our G fc,z -constrained and G p,9 /G fe,i -constrained multi-degree 
reduction algorithms. In each case, we give the least squares error E := y /and the maximum error 

:= max || P n (t) - R m (t)\\ « max ||P„(f) - R m {t )||, 
t£D N te[o,i] 

where Dn '■= {0, l/N, 2 /TV, ..., 1} with N := 500. 

In our experiments, we consider the “natural” choices for the values of parameters a, j3, i.e., (a,/3) £ 
{(0, 0), (|, |) , (—|) , (|, — |) , (— |, — |)}, and set the lower bounds d 0 , di of to 10“ 4 (see (14.181) 1. 

Taking into account the different types of continuity constraints, we compare the following cases: 

(i) G fc,z -constrained case (see m), which can be solved by using Theorem 14.21 

(ii) G p,9 /G fe,i -constrained case, solved by the algorithm described in > 15.21 

(iii) G fc,i -constrained case, solved by Algorithm 15.II 

Results of the experiments have been obtained on a computer with Intel Core i5-3337U 1.8GHz pro¬ 
cessor and 8GB of RAM, using 32-digit arithmetic. Maple™ 13 worksheet containing implementation of the 
algorithms and tests is available on the website http://www.ii.uni.wroc.pl/~pgo/GDegRed.mws. We use 
Maple™ f solve procedure, in the C p ’ q /G k ’ 1 case, to solve the system of linear equations, and QPSolve, 
NLPSolve procedures, to solve the quadratic and nonlinear programming problems, respectively. QPSolve 
uses the iterative active-set method, and for NLPSolve we select sqp method. Initial points for both proce¬ 
dures correspond to the values of continuity parameters in the C k ’ 1 case. 

Example 6.1. First, let us consider degree eleven Bezier curve which is an outline of the font “alpha” (for 
the control points, see [3, Example 6.1]). The results of multi-degree reduction are given in TableQ] Figs.flal 
and 111)1 illustrate two of the considered cases. One can see, that when it comes to minimizing E^ error, 
usually a good choice is a = /3 = — As expected, solution to the G k ’ 1 case is the most accurate, while 
C p,q /G k ’ 1 approach gives less precise results. C k ’ 1 conditions seem to be too restrictive, especially for k or l 
exceeding 2. 




Figure 1: Multi-degree reduction of degree eleven curve (blue solid line) to degree seven with C k ’ 1 (black dotted line), C p,q /G k ' 1 (green 
dash-dotted line) and G k ' 1 (red dashed line) continuity constraints; parameters: (a) a = (3 = — , p — q = 1, k = l = 2, and (b) 
a = j3 = —^,p=l,q = —, k = 3 and l = 1. 
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Parameters 



C fc,i solution 

CP’i/G k ’ 1 

solution 

G k ’ 1 solution 

m 

k 

l p 

9 

a 

P 

77 ( 0 ,/3) 
^2 

Eoo 

rj 2 

Eoo 


Eoo 

7 

2 

2 1 

1 

0 

0 

3.73 

5.97 

2.83 

5.27 

0.93 

2.26 





1 

2 

1 

2 

5.75 

5.83 

4.40 

5.14 

1.67 

1.91 





1 

2 

1 

2 

3.83 

7.53 

2.62 

6.42 

0.79 

2.91 





1 

2 

1 

2 

3.83 

7.69 

3.18 

5.18 

1.26 

2.19 





1 

2 

1 

2 

2.43 

6.10 

1.83 

5.40 

0.53 

2.54 

7 

3 

1 1 

- 

0 

0 

9.13 

16.24 

2.51 

5.11 

1.02 

2.41 





1 

2 

1 

2 

13.97 

16.72 

4.07 

4.95 

1.81 

2.07 





1 

2 

1 

2 

9.41 

19.51 

2.18 

6.38 

0.75 

3.11 





1 

2 

1 

2 

9.17 

18.79 

2.98 

4.91 

1.45 

1.99 





1 

2 

1 

2 

6.00 

15.82 

1.56 

5.25 

0.59 

2.69 


Table 1: Least squares error and maximum error in multi-degree reduction of degree eleven Bezier “alpha” curve. 


Example 6.2. Let us apply the algorithms to degree thirteen Bezier “heart” curve (for the control points, 
see [l2|, Appendix B]) and consider the case of k = l = 2. The results of experiments are given in Table [21 
Notice that the case of a = /3 = 0 was also considered in 12, §5.2] and 0, Example 4], As in 17|, we can 


clearly see that the solution to the G 2,2 case, in this paper obtained by Algorithm 15.II is more accurate than 
the result given by the approach proposed in [ 12 ], which leads to the G 1,1 /G 2,2 case (the same as for the 
algorithm discussed in *15.211 . As our approach considers different weight functions, it can be seen that the 
best choice to minimize Eoo error is a = /3 = — Fig. [2] presents a = /3 = 0 case. 

Now, we focus on the running times. We have implemented G 1,1 , G 1 , 1 /G 2,2 and G 1 ,:L /G 3,3 -constrained 
methods from 121, G 1,1 and G 1 ,:L /G 2 ’ 2 -constrained methods from Q as well as G 2,1 and G 2,2 -constrained 
methods from |l7}. The methods of Rababah and Mann and of Lu solve the same problem and give the 
same results as our G p, 9 /G fc ’-constrained method (see 1)5.21) . In Table [21 we compare the running times 


of these algorithms. Clearly, our approach is the fastest one. For the comparison of the G fe,z -constrained 
algorithms, see Table [4] Notice that, in some cases, our G fc,i -constrained approach is slig htly faster than 
the methods from 0 - We use Maple™ f solve procedure to solve the cubic equation fl7T ('2311 associated 
with the G 2,1 -constrained case. The implementation of G 2,2 -constrained method from 17] requires the 
unconstrained nonlinear programming solver. According to our experiments, the nonlinear simplex method 
(NLPSolve command with option method = nonlinearsimplex and the initial point A = r/ = 1) is the fastest 
solver available in Maple™ 13. Therefore, we use this solver for the purpose of the comparison. It is worth 
mentioning that Zhou et al. have omitted the constraints (14.1811 . Consequently, in some rare cases, the 
resulting curve may not preserve the original tan gen t directions at the endpoints. To avoid this issue, one 
can implement the improvements proposed by Lu [ 8 | . 


Parameters 

C 2,2 solution 

G 1 , 1 /G 2 ’ 2 

solution 

G 2,2 solution 

m a 

P 

771 ( 0 =, /3) 
rj 2 

Eoo 

E { 2 a ’ 0) 

Eoo 

TTiCo,/^) 

J ~ j 2 

Eoo 

8 0 

0 

1.52 

2.52 

0.64 

1.12 

0.36 

0.71 

1 

2 

1 

2 

2.37 

2.39 

1.05 

1.00 

0.62 

0.53 

1 

2 

1 

2 

1.58 

3.34 

0.64 

1.55 

0.42 

0.94 

1 

2 

1 

2 

1.52 

3.48 

0.74 

1.31 

0.37 

0.89 

1 

2 

1 

2 

0.98 

2.64 

0.39 

1.24 

0.21 

0.90 

Table 2: Least 

squares error 

and maximum 

error in 

multi-degree reduction of degree thirteen Bezier ‘ 

‘heart” curve. 
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Figure 2: Multi-degree reduction of degree thirteen curve (blue solid line) to degree eight with G 2 ’ 2 (black dotted line), C' 1,1 /G 2,2 
(green dash-dotted line) and G 2 ’ 2 (red dashed line) continuity constraints; parameters: a — (3 — 0. 



Parameters 



Running times [ms] 


m 

k 

/ 

P 

Q 

Our C p ’ q /G k ’ 1 

method Rababah and Mann [17] 

Lu [6] 

8 

1 

1 

- 

- 

32 

62 

63 

10 

1 

1 

- 

- 

31 

94 

63 

12 

1 

1 

- 

- 

47 

94 

62 

8 

2 

2 

1 

1 

16 

63 

47 

10 

2 

2 

1 

1 

31 

63 

62 

12 

2 

2 

1 

1 

31 

78 

78 

8 

3 

3 

1 

1 

15 

78 

— 

10 

3 

3 

1 

1 

31 

109 

— 

12 

3 

3 

1 

1 

63 

no 

— 


Table 3: Running times of C p ’ 9 /G fc,i -constrained multi-degree reduction of degree thirteen Bezier “heart” curve; parameters: a = /3 = 0. 


Parameters Running times [ms] 

m k l Our G k ’ 1 method Zhou et al. [17] 

8 2 1 92 108 

10 2 1 121 137 

12 2 1 168 166 

8 2 2 153 204 

10 2 2 298 248 

12 2 2 290 292 

Table 4: Running times of G fe,i -constrained multi-degree reduction of degree thirteen Bezier “heart” curve; parameters: a — [3 = 0. 
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7. Conclusions 


In this paper, we propose efficient methods of solving the problems of G—constrained and C p,q /G k ' 1 - 
constrained multi-degree reduction of Bezier curves with respect to the least squares norm. We give two-phase 
algorithms of solving these problems. 

The first phase of the algorithms consists in finding values of the geometric continuity parameters to mini¬ 
mize the error m■ In the case of G fc —constrained problem, we solve the quadratic or nonlinear programming 
problem to obtain these values. For G p,9 /G fc —constrained case, we use some simplifying assumptions, i.e., we 
impose constraints of G-continuity at t = 0 when k > 1, and at t = 1 when l > 1. Therefore, by fixing some 
of the parameters, this approach leads to the system of linear equations ([4. 91) . Assuming that — 1 < k,l < 3 
and including the hybrid cases, there are 37 continuity cases which require computation of the continuity 
parameters. Those variants of the problem differ, and we have not proven that in each case a unique solution 
exists. 

During the second phase, which is the same for both approaches, we use the properties of constrained dual 
Bernstein basis polynomials to compute control points of the multi-degree reduced curve. The complexity of 
this phase is O(mn), where n and m are the degrees of the input and output curves, respectively. This is 
significantly less than complexity of other algorithms. Moreover, our approach avoids matrix inversion. 

As expected, solution to G—constrained problem is the most accurate, while the one obtained by 
C p ’ 9 /G fc ’—constrained multi-degree reduction is less precise. C k ’ 1 conditions tend to be too restrictive, es¬ 
pecially for k or l exceeding 2. Comparison of running times of our C p ’ 9 /G fc ’—constrained approach with 
analogous methods from [6:, |l2| shows advantage of our algorithm in practice. Furthermore, the experiments 
show that our G fc —constrained approach is comparable to the methods of a even slightly faster in some 
cases. 
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